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START: 



Set Up 
Parallelization 



Initialize 
Constants and 
Sine/Cosine Tables 



Input 1 & 
Initialize the List of 
Fourier Indices (hkl) 



Get Biggest Sphere 

Wiihoul Overlap 
of Symmetry Mates 



□ Input 2 



Confirm 
Cipfciilation Mode 



■ Initialize 
Ta|B|es of Spherical 
BeRsel Functions and 
Bessel Function O's 



Determine Limits 
on 1 & m indices and 
on n Indices at each 1 



For modes 5 thru 7: 
Read expected SHSB 
coefficient & 
deviation values 



For Use with a Radial Correction or 
with Modes 5 thru 7: 



Mode$l&2{contd) 



Get 

Fractionalization 
Matrix 



Partition the List 
of Fourier Indices 
Between 
Proce?>sor(s) 



Convert Each Fourier 
Index (hkl) to Polar 
Coordinates &. Get 
Its Bcssci Argument 



Get the FT of a 
Ciystallinc Unit Cell 
Filled with Symmetry 

Related Spheres 



Calculate a Radial 
Correction Factor 
for Each "n" Index- 



Calculation Mode-Specific Routines: 

Modes 1 & 2 
(Unpliased Diffraction Amplitudes to 
Phased FT of SHSB-modeled Unit Cell) 



Modc-SpccifiC 
Input 



Get Input File # for 
Calculated FT, if 
there is a Prior Model 



Initialise Correlation 
Exponentiation 
Factor 





Cycle 




For Each Choice (mstop) of 


Stopping Value for the SHSB 




m index 






Update Correlation 






Exponentiation 






Factor 





If there is a Prior Model 
(On cycles after the 1st cycle 
of the 1st run of the program): 



On Cycles After the 
iRt value of mstop 

Update the File # for 
the Calculated FF 



For Each input SILSB M odel 

Convert the File # 
& the Model # 
to a File Name 



Read Input 
FT & SHSB coefs 
of the Prior Model 

"lf"there Is No Prior Model: 



Model the Crystal as 
a Crystalline Unit Cell 
Filled with Symmetry 
Related Spheres 



JUL 



Get a Starting Value 
(msus) for the SHSB 
m Index for the Next 
Packet of m Values 



Update the Model 
by Adjusting the 
Contributions of 
each SHSB function 
for each Imn index in 
the m-pucket from 
'msus' to *mstop' 



Fkw Chart for tbo Main Driv^ Options toconi^HitDatheFrof aSHSBModelof 



Figure T A schematic example: Two choices for filling the same 
portion of a crystal unit cell from an orthorhombic Spacegroup, 
Although the spheres on the right are smaller than those on the left, 
for some crystals, the local maximum in the packing on the right 
wold be the packing of maximal consistency with the 
crystaliographic data. 



Initialize 
Fractlonalization 
Matrix 



Initialize the Equal 
Partitioning of the 
Fourier (hkl) Index 
between Processors 



On Ist Cycle of 1st Run ; 

Prescale Observed 
Diffraalon to that 
of a Unit Cell of 
Spheres 



^^Sefine the First 

ex. 

StjSB Index Triplet 
(ti|in)for which to 
CpK^sider Model IFj's 




!i^|tiallze for Index- 
:i^-lndex Update 
bf^Origin-Centered 
SplSB Basis Function 




Mibdes 4 & 5 only: 
l|)|lallze Buffers for 
CUffl^ulative Update of 
FoiJifjer Representation 



Initialize Pointers to 

Stored Fourier 
Representations of 
Model and of Basis 



Mode 3 only; 

Get File Name fronn 
File #& Open It to Let 
SHSB Coefs^be Read 



1 For each "nn" Index 
1 (0 to maximum "m") 




For Each hkl in this 
Partition: 

Update "pf)" Recursion 
Formula for Fourier 
Representation of the 
Origin-Centered SHSB 





For each 1" index 
(present "m" to maximum "1") 




For Each hkl in this 
Partition: 

Update "l;m" Recursion 
Formula for Fourier 
Representation of the 
Origin-Centered SHSB 





For each "n" Index 
(1 to maximum "n" for each "1") 




For Each hkl In this 
Partition; 

Update "n" Recursion 
Formula for Fourier 
Representation of the 
Origin-Centered SHSB 





Depending on Mode: 
Choose the # of Passes 
and # of Presumed 

Phase Angles Needed 
for the SHSBcoef. with 

this SHSB index (Imn) 



Set the Presumed 
Amplitude of the 
Origin-Centered 
SHSB Basis Function 



FIRST PASS: 

Initialize Registers: 
Overall Comparison 
of Correlation Coef 
& Other Statistics 

Renitialize Pointers to 
Storage Sites for Fourier 

Representations of the 
Full-Unit-Cell SHSB Basis 

Parallel Processor Version: 

Set # of calculations to: ' 
(# of presumed values of ^ 
SHSB coef.'s phase) , 

X ! 

(# of stored accumulated ' 

SHSB models for trial ' 

combination v^ith this , 
new SHSB component) 



Given: # of processors 
# of hkl partitions 
# of calculations 
Get:# of required 
rounds of trial 
combinations 

For each round of trial 
combination on this processor 

Single Processor Version: 
(Outer Loop) 

I For each presumed value 
of the SHSB coel's phase 



Initialize Registers: 
Angular Comparison 
of Correlation Coef. 
& Other Statistics 



Via J- 



